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q ; A new model for calculating the structure of bound states of interacting 

particles is considered. The model takes into account the noncommuta- 
tivity of the space and impulse operators plus the correlation equations 
for the indeterminacy of these quantities. The efficiency of the model is 

oo 1 demonstrated by specific calculations for some lightest nuclei. 
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>'. 1. Introduction. The formulation of the model. 
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In our recent work (Tj we have proposed a new quantum mechanical model for 
interacting bodies. The idea [2J- that the coordinate and impulse operators for 
different particles may be not commutative make up a basis of the NOCE model 
^ ; (the NOCE model means the noncommutativity of the operators and the correlation 
equations). Within the framework of the NOCE model, we have examined [lj in 
detail the ground and some excited states of Hydrogen-like (if-like) atoms. In the 
present paper the 2 H — D, 3 H and 3 'He lightest nuclei are studied by a given model. 
In the case of A identical particle the evident generalization of the equations 
^ ! (5)- (8) from the paper [lj yields the following relations 

_ _ i 

Mn = [y k ,P y i] = [z k ,f l } = l fA 1 l = J. (i) 



k^l; 

pi = -ih/3 ■ \/ b k,l = 1,2,..., A; 

where there are only two noncommutativity parameters (3 and j3 present. Assuming 
the commutator of the operator of the particle coordinates and the operator of the 
total impulse P to be equal i% one obtains from the above mentioned equation a 
simple connection between (3 and (3 

l = /? + (A-l)fl,. (2) 
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For this reason, only one additional equation is needed in order to find the noncom- 
mutativity parameters (3 and j3 . We write this equation by examining the process 
of the measurement of the coordinate of k-th particle with maximum accuracy. The 
reasoning similar to that employed in [T] leads us to the equation connecting the non- 
commutativity parameter (3 and the matrix element (ME) of the force f a = (\F\), 
exactly, 

Po = -=~lo • f ; £ = mc 2 . (3) 

£ z 

Here m is the nucleon mass, and 7 G is one of the correlation factors specified by the 
equations 
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P, k = l; 
■ P , k y£ /; 



(4) 



Ax k Apf = Ay k Apf = Az k Apf = 

/c, / = 1 , 2 , . . . , A. 
Eqs.(2) and (3) yield one more useful relation : 

P = 1- -^-(A-lhofo. (5) 

£ L 

In the NOCE model, we have the analogy Schrddinger equation (SE) 



h 2 A V72 1 



7 E VJ + .E, V(Tij) 



2m i=l - i>j=l 



E 



(6) 



m' = m/ l-%-(^-l)7 /, lL> 



Specific solutions to these equations can be found by the method of successive iter- 
ations, where at the 1-st step the conventional SE with the masses m = m, which 
the particles have in the absence of the interaction {(3 = 1), has to be solved. After 
that, on finding the wave function ijj, one can calculate ME of the force and 
the first value of the commutation parameter (3 distinct from unity. At the 2-nd 
step, SE is solved with the modified particle masses m = m/P 2 . On finding the new 
ijj, we calculate the quantity (\F\) and compare it with the one obtained at the 1st 
step. Then, we proceed with the iterations until the values of the matrix element 
of the force (\F\) obtained at subsequent steps will be virtually indistinguishable. 
It is clear that, before starting the above iteration process, we should specify the 
numerical value for the correlation factor 7 entering Eqs.(6). To calculate its, one 
can employ specific parameters of a given system based on reliable experimental 
data. The way to practically implement this will be described in detail hereinafter 
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2. Deuteron and the variational Ritz principle. 



Now we proceed to consideration of nuclear systems taking as an example the 
simplest one, the deuteron. As is known, in contrast to the atom theory, in the 
nuclear theory one has to deal with a serious problem of choosing the nucleon-nucleon 
(NN) forces. In view of this circumstance the accuracy of theoretical estimates for 
a nucleus is considerably lower than that of similar calculations in atomic physics. 
For this reason, taking into account the calculations of i7-atoms given in [Ij, such a 
small correction in relation to the results of conventional SE may evidence, at first 
sight, the inefficiency of application of the NOCE model to nuclei. However, if it 
is remembered that the matrix element (ME) of the force {\F\) for nucleon systems 
should be several orders greater than similar values for (\F\) in atoms, the contrary 
anxiety arises, meaning that the application of the NOCE model to nuclear physics 
may result in definitely unrealistic theoretical estimates. 

Apparently, to clarify these questions, one needs to perform specific numerical 
calculations. First, we can take as an example the deuteron with a simple central 
potential in the Yukawa form 

p -r/a 

V(r) = -V —. (7) 
r/a 

The bound states of the deuteron are described by the wave function ^(r) = 
-Xni(r)Yi n (6 , if) depending on the vector of the relative motion r = Y\ — Y2- Let us 
restrict ourselves to considering the ground s-state (/ = 0), then, the radial wave 
function Xn=i,i=o{ r ) = x{ r ) nas to satisfy the simple equation 



h 2 d 



2 



+ V(r) - E 



2 // dr 

where the reduced mass \i modified in the NOCE model equals 



XW = 0, (8) 



/ llb-l 1 1 I/O IIC 9 

II = — ; r = o ! £ = TTIC . 9 

mi + m 2 2(l-^ 7o (|F|); 2 



To write the equation (9), we use the relation 

TYlproton ^neutron — Wl i (1^) 

wich are natural for the case of identical particles. Satisfactory solutions to the 
equation (8) can be found by means of the variational Ritz principle. For the 
Yukawa-type potentials, good results can be obtained by using the function 

X (r) = 2r/ 3/2 re^ r (11) 
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as the variational one. Here the only variational parameter rj has to be found from 
the condition that the deuteron energy 

must take extremum. In what follows, it is convenient, following Ref. |4J , to introduce 
the notation 

J? 

Then the above mentioned condition dE/dr] = yields the cubic equation with 
respect to the quantity p 



K=-^-V , p=2ar 1 . (13) 



p 3 + (3 - K) p 2 + 3(1 - K) p + 1 = . (14) 

As is evidenced by the immediate calculation, the optimal value is associated with 
the root p , which can be expressed via the quantity K introduced in (13) as 



K 

p = — -l + 2> 



K( K\ 7T 1 

— 11 + — ) • cos < — — - arccos 



ts K 2 K 3 
^ 6 27 

(f + f ) 3/2 



(15) 



3 iT l 3 \ 3) w °l3 3 
The matrix element of the force 

F(r) = ~V{r) = -aV (\ + —) exp(--) (16) 
dr \r z arj a 

on the functions (11) is equal 

/„= (x(r)\F(r)\ X (r)) =&\a V f \ = ^. p "jl±SL (17) 

Substituting the value of the root p Q from (61) instead of p into the relation (17) 
and taking into account that the quantity 

is, in turn, expressed via f a: we obtain some rather complicated equation for ME 
of the force f Q . The equation found in this way plays, apparently, the same role as 
the equation (33) from (Tj does in the case of i7-like atoms, namely, solving this 
equation enables us to avoid the complicated (in computational respect) procedure of 
successive iterations. The necessary solutions have been found in the graphical way. 
In doing so, from the set of all the solutions found (i.e., the points of intersections of 
the straight line y = x = f with the plot of the function y = (fi2{x) representing 
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the right-hand side of the equation (17)) we selected the one associated with the 
minimum value of f > 0. 

In the numerical calculations performed here, we used the parameters of the 
potential (7) given in [5J, V = 20.7 MeV, a = 2.43 fm. This potential was adjusted 
to fit the experimental data on the scattering length and effective radius of the 
triplet state in the "proton + neutron" system. The nucleon mass was equal e = 
mc 2 = 931.441 MeV. The calculations were carried out by using the various allowed 
([0 < 7 < 1]) values of the correlation factor j . The results obtained are given 
in the Table 1. In this table, the ground state energy of the deuteron E = Ed, 
ME of the force / , the optimum value of the root p = 2a t/ , and the values of the 
quantities 



defining the commutation relations (1) are presented as dependent on the quantity 



As is seen from the table, the binding energy of the deuteron \E^\ (the exper- 
imental value \E^'\ = 2.22457 MeV) grows with 7 from \E D \ = 2.2019 MeV at 
7 = 0.00001 up to \E D \ = 2.3826 at 7 = 1. This means that in the NOCE model 
under consideration, for the potential of TViV-forces chosen, the maximum possible 
increase of \Eu\ constitutes ~ 8% as compared to the calculation of \Ed\ based on 
the conventional SE. With increasing the 7 G parameter from 0.00001 to 1 the quan- 
tity (3 diminishes insignificantly (by less than 1%) while the quantity (3 grows from 
« 0.906 • 10~ 7 at 7 = 0.00001 to (3 « 0.01 at 7 = 1. 

The experimental binding energy of the deuteron corresponds to the values 7 = 
0.13523, f = 20.0812MeV//m, and (3 = 0.998 765, accompanied by the increase in 
the nucleon mass by the value of Am = {-h — 1) m = 0.00247 516 m. 

3. Mirror nuclei 3 H and 3 He. 

In this section, we consider the results of the calculations of the ground states 
for the Tritium and Helium-3 nuclei. 

In the case of a central iViV-potential V(ry), the force F entering the definition 
for the quantity f is equal to the sum of derivatives of V(rij) with respect to the 
relative distance ty,-, i.e. 



l-2x = /3, 2 



o — 




(19) 




r 



(20) 
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In specific calculations of light nuclei, one frequently restricts consideration to the 
case of the central exchange potential 

V { r ij) = ~ E V 2 S+l,2T+l(nj) As+l,2T+l(u) , (21) 
S,T=0,1 

where the known projection operators A^+i^t+i do "cut" the relevant 
spin-isospin states of the interacting nucleon (ij)-pair from the wave func- 
tion ¥(1, 2, A). The radial dependence of the components of iViV-potential 
V2S+i,2T+i{rij) can be presented, without loss of generality, in the form of a series in 
Gaussian terms 

V-pot / f^.\ 

V 2 S+l,2T+l(r l j) = E ^25+l,2T+l • ex P — | • (22) 
v=l V AW 

Now we consider the problem of the bound states for the nuclei 3 H and 3 if e. In the 
NOCE model, the conventional three-body SE 

#¥(i |2| 3) = # = Ef-^W, 2 + -^-+ E n|ri-r 2 |), (23) 

i=i V 2m j |ri - r 2 | 1>3 = X 

with the modified nucleon mass 

m=- 2 , /J = l-_ 7o / . (24) 

is valid. The charge entering (23) is Z = for 3 H nucleus (1-st and 2-nd particles 
represent neutrons), and Z = 1 for the nucleus 3 He ( 1-st and 2-nd particles are 
protons), respectively. 

In specific investigations of the structure of light nuclei, we use as the iViV- 
potential of the type (21)-(22) the two-Gaussian Volkov potential |6J with the 
following parameters (I-st version), 1-st component {y = 1) vji = V-}^ ^ = 
144.86 MeV , V^ 1] = v}" =1] = -28.972 MeV, ^1=0.82 fm; 2-nd component {y = 2) 
V^ =2] = V}r 2] = -83.34 MeV, V^ 2] = V^ =2] = 16.668 MeV, /i 2 =1.60 fm. 

In view of modern state of computer technology, it is convenient to seek the 
solution of three-body SE by using the variational Ritz principle with the wave 
function of the system being expanded in a series of the known basis functions. These 
latter may contain one or another set of variational parameters, which, actually, 
determines the flexibility of the basis chosen. In the case of the lightest nuclei, 
satisfactory results were obtained by using the simplest Gaussian-type functions. 
The efficiency of using the Gaussian basis of functions in few-body problems has 
been first realized by the theorists of the Moscow State university in 1973-1975 (see, 
for instance, [7j - (ID])- Really, even the first calculations performed within the 
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framework of the stochastic variational method (SVM) gave the results of highest 
accuracy possible at that time. In view of this, some authors sometimes refer to 
these calculations as the "high precision" ones. Nowadays, the SVM method is 
widely and successfully employed in theoretical studies of various quantum systems 
consisting of small number of particles [TT] - [13j. 

That is why, it is reasonable to employ here SVM to calculate the lightest nuclei 
under consideration. In doing so, the total wave function *H and 3 He is presented, 
as usually, in the form of a product of the space and spin-isospin function 

^(l,2,3) = ^(n,r 2 ,r 3 )-x(l,2,3). (25) 

According to the experimental situation, the function \ represents the determinant 



X(l,2,3)= ' 



xtvK-0 xt,t( 2 ) xt,t( 3 ) 

Xj,^ 1 ) X|,t( 2 ) X|,t(3) 
Xut 1 ) Xu( 2 ) Xu( 3 ) 



(26) 



composed by the one-particle spin-isospin functions 

where the lower subscript m a of the function (27) takes two values denoted in (26) by 
an arrow "up" f or "down" |. Exactly same notation is used for the other quantum 
number m T . Resolving the determinant makes it evident that the function 2, 3) 
is antisymmetric with respect to permutations of spin-isospin coordinates of 1-st 
and 2-nd particles. For this reason, by virtue of Pauli principle, the ^-function must 
be symmetric with respect to permutations of space coordinates ri and r 2 , which 
means, in turn, that the expansion of ^-function should be carried out over the 
properly symmetrized Gaussian terms, i.e., 

</>L=o(ri, r 2 , r 3 ) = E Cj • Uj(r 1: r 2 , r 3 ) , 

Uj(r u r 2 , r 3 ) = A \j) = A exp{-a J 12 r\ 2 - a{ 3 rf 3 - a J 23 r| 3 } = (28) 

2 

= E exp {-a{ 2 (sym) ■ r\ 2 - a{ 3 (sym) ■ r 2 u - a J 23 (sym) ■ r| 3 } . 

sym=l 

The action of the symmetrization operator A on the usual Gaussian term \j) gives 
rise to a sum consisting, in this specific case, of two summands only with param- 
eters a J kl being replaced by a J k i(sym). The result of the immediate examination of 
permutations associated with the operator A is given in the Table 2. 
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As mentioned above, in numerical calculations we used the central exchange 
Volkov potential, which contains the dependence on the spin-isospin coordinates in 
the projection operators only. This circumstance makes it possible to easily calculate 
the matrix elements (ME) of the operator of NN -forces (21) with the spin-isospin 
functions. In particular, for Volkov potential in the case of the nuclei 3 H and 3 He 
under consideration, we obtain 



(X(l,2,3) 



i>]=\ 



X(l,2,3)}= (29) 



= V 13 (r 12 ) + \ [3V 3 i (r 13 ) + V n {r 13 ) + 3^23) + V n {r 23 )} . 

Now, in order to write the system of linear equations for the expansion coefficients of 
the wave function Cj : we calculate the matrix elements of all the operators entering 
the definition for the Hamiltonian (23) on the basis function Uj(r\, r 2 , r 3 ). Earlier, 
the analytical calculations of this kind have been performed in the number of works 
(some of them are cited above). In view of this, we cite here only selected results 
sticking to the notation of the paper by N.N.Kolesnikov [10J. Let us begin with the 
overlap integral for the basis functions : 

2 2 / / 7T \ 3 

(uy\uj) = £ E (j,zym\j,sym); (j , zym\j , sym) = [-==) , (30) 

with the notation 

\h s V m ) = ex P {-a{ 2 (sym) ■ r\ 2 - a{ 3 {sym) ■ r 2 3 - a J 23 {sym) ■ r| 3 } ; 



(31) 



,j\ = a 12 a 13 + a u a 23 + a 13 a 23 ; a k i = a J kl (sym) + a 3 kl (zym) 



It is clear that ME for any other operator on the functions (28) represents a sum of 
the type (30), i.e, the analytical calculations reduce, actually, to finding the partial 
ME with the functions |j, sym). For the operator of the kinetic energy, we obtain 



(j',zym\f\j,sym) = ^ «|«) . 



3 3' 

E E a J kl (zym)a J kl (sym) 
fc=ij',i=i 



U l'k + U lk 



J J\ 



(32) 



with the notation [10J 



D 



[3 j] — d 



Ik — da M 



Dti J \ D 



D 



[3 3] _ 



12 



3 3] _ 



r ■' -l 

«13 + «23, ^13 = «12 + «23, ^23^ = «12 + «13 



0; 



(33) 



3 j] _ 
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The calculation of ME of the operator of potential energy is associated with the 
expression 

\ 3 



/ 



Vpot 

(j ,zym\V(rki)\j,sym) = J2V [ 



7T 



(l',l^k), (34) 



Vpot 



where 

V(r kl ) = E^ M -exp(-^f| (35) 

v=l \ H v J 

is one of the components of the iViV-potential (21) remained in (29) after averaging 
over spin-isospin functions x(l,2, 3). When considering the nucleus 3 He, we have 
to calculate also the Coulomb ME 



(j\zym 



j, sym) = ^= (j, zym\j, sym) 



(36) 



Now, we have all the relations needed to write the system of equations 

H ■ X = XB ■ X , 



(37) 



which is necessary to find the energy spectrum of the nucleus {Ei} = A and the rele- 
vant coefficients in the expansion {Cj} = X of the space wave function ^>(ri, T2, r^j. 
According to the definition of the matrix elements B — \\B^j\\ and H — \\HsA\ in 
Eq.(37), we obtain 



j j 1 



U f \Uj 



= E E (j \zym\j,sym) , 

zym=l sym— I 



(38) 



2 2, 

E E (j , zym\H = T + U + C/conZ.li, s ?/ m ) • 

zym=l sym=l 



Thus, the energy of three nuclei is sought as a result of solving, in the first place, 
the generalized problem for eigenvalues, and, in the second place, the problem of 
finding the optimum values for the variational parameters a 3 k i determining the basis 
functions (28). In so doing, the most difficult procedure herewith is the optimiza- 
tion. In performing numerical calculations, we employed the known mathematical 
library IMSL, specifically, we used the subroutine DGVCSP [Hj in the generalized 
problem for eigenvalues, and the subroutine DBCONF [T^j for the optimization 
procedure. 

It should be kept in mind that the above mentioned calculations are carried 
out only after the value of the nucleon mass m = m/fi 2 is found. The latter is 
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determined by the magnitude of ME of the force f Q (see the relation (24)). In this 



f =(*(!, 2, 3) 



dV{n 



3 

^ ilr 

i>j = l 



tf(l,2,3)) = (39) 



= wxmwm X c i c i X {v * ■ G ' /jl < 1 - 2 )+ 

+i(3^r ] + Wr ] ) ■ g/ j1 (i,3) + &3v& ] + wr 1 ) • ^(2,3)} , 

with the partial ME being 



— 142 



j, sj/m) = (40) 




The quantities ^ and d\{^ entering this equation have been already specified 
by the relations (31) and (33), respectively. On calculating ME of the force f in 
this way, we evaluate the mass m by Eq.(24) and proceed to solving SE with the 
modified mass m , i.e., we carry out the process of successive iterations described 
above (Section I). In so doing, we may put f o = (m = m) at the first step, then, 
the 1-st iteration is identical to solving the conventional SE. 

Consider now the calculations of the nuclei 3 H and 3 He carried out for Volkov 
AW-potential JB] with the model parameter 7 G = 0,6. In this calculation, the 15 
functions (28) were employed. As is seen from the definition (28), each of these 
latter contains three independent variational parameters oj 2 5 a i3 an d 0^3. Some of 
the results of the calculations performed here are given in the Table 3. The values 
of the energy and ME of the force f are given in the units of MeV and MeV/fm, 
respectively. Since the iViV-potential used in calculations was adjusted to fit the 
basic properties of light nuclei, the calculated binding energy EQ- mech - turned out to 
be close to the experimental value. As one would expect, the NOCE model overbinds 
the nuclei 3 H and 3 He to some extent (see Table 3). The energy gain constitutes 
~ 4% as compared to E^' m , and the change in the nucleon mass is ~ 0.9%. 

It is clear that these numbers directly depend on a specific choice of the AW- 
potential, as well as on the magnitude of the parameter For this reason, the 
theoretical results presented here are qualitative in nature, in contrast to atomic 
calculations [1J. However, they are sufficient to make estimates of the efficiency of 
the NOCE model under consideration. In others words, the calculation performed 
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made possible the estimation of the order of the expected corrections to the basic 
nuclear properties arising due to the allowance for the noncommutativity of the 
coordinate and impulse operators of the interacting particles. 
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Table 1: The properties of the ground state of the deuteron within the NOCE model. 
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Table 2: Symmetrized variational parameters {a J kl (sym)} for the nuclei 3 H and 3 He. 
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Table 3: The properties of the ground states of the 3 H and 3 He in the case 7 C = 0.6 and the Volkov 
potential. 
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